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Abstract — Matrix theory and its applications make wide use of the eigenprojections of square 
matrices. The paper demonstrated that the eigenprojection of a matrix A can be calculated 
with the use of any annihilating polynomial for A", where u > indA. This enables one to 
establish the components and the minimum polynomial of A, as well as the Drazin inverse A^. 

1. INTRODUCTION 
1.1. On Applications of Drazin Inverse, Eigenprojection, and Matrix Components 

The Drazin inverse, eigenprojections, and components of square matrices, as well as their mini- 
mum polynomials, are widely used in mathematics, mechanics, engineering, economics, and other 
sciences. Let us mention some of these applications. Dynamics of multibody systems is described 
by systems of differential-algebraic equations. Their linearization results in "implicit" singular sys- 
tems of ordinary differential equations. Their analysis relies on determining the Drazin inverse 
for the matrices of coefficients [1]. This approach was used in [2] to analyze stabilization of the 
model of truck motion. The eigenprojections of a matrix that correspond to its eigenvalues are 
its main components. From the computational point of view, determination of the eigenprojection 
corresponding to the eigenvalue amounts to determining its Drazin inverse and can be used to 
determine the Moore-Penrose inverse. The matrix components and the functions of matrices de- 
termined with their use are employed in the dynamics of rigid bodies [3], motion stability analysis, 
and integration of systems of linear differential equations [4,5]. Applications of the eigenprojection 
and Drazin inverse in the theory of Markov chains [6], graph theory [7], cryptography [8], and 
econometrics [9] (see also [10]) deserve mentioning as well. In all of these applications, one needs 
to determine eigenprojections, Drazin inverse, and matrix components. The present paper shows 
that for this purpose one may use any annihilating polynomial for j4", where j4 is a given matrix 
and u > ind^. 

1.2. Basic Definitions 

Let A G C"X" be a square matrix whose range and null space are denoted by 'R-{A) and 7V(j4), 
respectively. Let r be the total multiplicity of the nonzero eigenvalues of ^, v = n — r he the 
multiplicity of as the eigenvalue of A, and v = indA be the index of A, that is, the smallest 
k G {0,1,...} for which rank(A'^+^) = rank(A'^) [11,12]. The index is if and only if A is 
nonsingular. The index of a singular matrix is the index of its eigenvalue 0, that is, the multiplicity 
of zero as the root of its minimum polynomial, or the size of the greatest block with zero diagonal 
in its Jordan form. 

^ This work was supported by the Russian Foundation for Basic Research, project no. 02-01-00614-a. 
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By the eigenprojection^ of a matrix A corresponding to the eigenvalue [14] or, for brevity, just 
by the eigenprojection of A [15] is meant the projection, that is, the idempotent matrix, Z such 
that 

7^(Z)=A^(A^) and M{Z) = n{A''). 

Z can be said to be the projection on M{A^) along Tl{A^). The eigenprojection is unique because 
an idempotent matrix is uniquely determined by its range and null space (see, for example p. 50 
in [11]). 

The eigenproj actions are used in the theory of generalized inverses and numerous applications of 
linear algebra, as well as for decomposing matrices into components^ which allows one to determine 
the values of f{A) for a wide class of functions / : C ^ C (sec [3-5]). In this connection, we also 
note that for any finite homogeneous Markov chain with the transition probability matrix P the 
limiting matrix of mean probabilities 

1 k-l 

p°° = lim - y P* 

k^oo k 

is the eigenprojection of the matrix I — P [6, 14]. Another example: for any weighted digraph F 
with the Laplacian matrix L, the matrix P of maximum converging forests is the eigenprojection of 
the matrix L [7] . The present paper is devoted to the methods of determining the eigenprojection 
of a matrix, as well as its components and the minimum polynomial. 



1.3. Preliminaries 

From the fact that is the order of the greatest Jordan block with zero diagonal in the Jordan 
form of A, it follows that 

u<v. (1) 

For the powers of a Jordan block J with zero diagonal we get 

rank {J^) = max{Z — k, 0}, 

where I is the order of J. Therefore, 

rank (A*) > rank(A*+^) for < i < u and rank(A') =r for i > v and (2) 

ind {A'') = mm{teZ\kt>i^} = \u/k'\ , (3) 

where \x\ is the smallest integer not smaller than x. 

Since 7?.(X) and M{X) are subspaces of C" for any matrix X G C"^" and since 

dmiUiX) =n- dimJ\f{X) = rankX, 

we obtain by using (2) that 

Af{A) c ... cMiA") =M{A''+'^) = (4) 

n{A) D . . . D 7^(A^) = n{A''+^) = ..., (5) 

where the strict inclusions refer to the case oi v > 1. These relations and (3) give rise to the 
following lemma. 

Lemma 1. The eigenproj ections of the matrices A, A"^, . . . A'^ , . . . are the same. 
^ The eigenprojection is also called the principal idempotent (see, for example, [13]). 
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1.4- Plan of the Paper 

The paper makes use of the following approach. The eigenprojection of a matrix with index 
or 1 is the easiest to determine. Therefore, knowing an upper bound u for the index of the 
matrix A, one can seek by virtue of (3) the eigenprojection for which has index 1 instead 
of the eigenprojection for A. Indeed, according to Lemma 1, the eigenprojections of A and 
coincide. In the next section, we use this approach and several well-known characterizations of the 
eigenprojection to obtain simple expressions for the eigenprojection of and, consequently, of A. 
In Sec. 3 we propose a method for determining the eigenprojection of A (and the Drazin inverse 
A^) by means of any nonzero annihilating polynomial for A^. Section 4 discusses calculation 
of the components and the minimum polynomial of a square matrix. Explicit expressions of the 
eigenprojection and the components of a matrix A in terms of its eigenvalues are given in Sec. 5. 



2. CHARACTERIZATIONS OF THE EIGENPROJECTION 

Below we present a number of conditions that are equivalent to the fact that an idempotent 
matrix Z G C"^" is the eigenprojection of a matrix A € C"^" having index v. This list refers also 
to the publications where these conditions are presented. 

(a) TZ{Z) = MiA") and = Mi{A*Y), where X* is the Hermitian adjoint of X [14]. 

(b) A'^Z = ZA" = and rank^l'' + rankZ = n [16, 17]. 

(c) AZ = ZA and A + aZ is nonsingular for all a 7^ [18]. 

(d) AZ = ZA, A + aZ is nonsingular for some a ^ 0, and AZ is nilpotent [18]. 

(e) AZ = ZA, AZ is nilpotent, &ndAU = I-Z = VA for some U,V e C"^" [19]. 

(f) ZC = CZ for any C commuting with A, AZ is nilpotent, and (det A = 0^ Z j^O) [20]. 

Now we pass to some more constructive characterizations of the eigenprojection. In the items 
to follow, u is any integer not exceeding u, and we use Lemma 1 according to which A and A" 
have the same eigenprojection. Also, the following notation is used: A^ is the Drazin inverse for 
A, (^4")* is the group inverse for A^, 

n 

V^(A) = ^a„_,V (6) 

j=v 

and 

n-l 

B{X) = adj(A/ -A) = Y, An-i-jX^ (7) 

j=o 

are, respectively, the characteristic polynomial for A and the adjoint matrix for XI — A; ipuW and 
'0„(A) are, respectively, the characteristic and minimum polynomials for A^; du{X) is their quotient: 

xPuW = duiX)M\); (8) 

BuW = adj(A/ — is the adjoint matrix for XI— A^; and C„(A) is the reduced adjoint matrix [4] 
for XI — A^, that is, the matrix polynomial 

which can also be obtained by dividing the matrix polynomial ^uiX)I by the binomial XI — A". 
Some constructive characterizations of the eigenprojection are as follows. 
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(g) Z = I- AA^ =1- [13, 15]. 

(h) Z = I-{I- An^^/an^vT [13]. 3 

(i) Z = Cu(0)/p„(0), where Pu{^) is the polynomial obeying the equality Apu(A) = V'u(A) [4]. 
(j) Z = lim AS„(A)/V'«(A), where A : Vu(A) ^ 0. This follows from (8), (9), and (i). 

A — ^0 

(k) Z = lim (/ + ryl")-\ u>v. This follows from (j) or Theorem 3.1 of [22]. 

\t\^oo 

(1) Z = X{Y* X)~^Y* , where X and Y are the matrices whose columns make up bases in J\f{A'^) 
and M{{A*Y)^ respectively [13, 14] {v can be replaced by u by virtue of (4) and (5)). 

In the following section, we propose another constructive characterization of the eigenprojection. 
We do not touch upon the issues of computational efficiency and just note that the corresponding 
algorithm does not require computation of limits, the reduced adjoint matrix, the matrix coefficients 
of the adjoint matrix, generalized inverses or bases of subspaccs. If an arbitrary annihilating 
polynomial for A^ or the eigenvalues of A are known, then the eigenprojection of A can be expressed 
as a scalar polynomial of A^. If the characteristic or minimum polynomial is used as the annihilating 
polynomial, then the resulting expression is equivalent to the expressions presented in (j) and (i), 
respectively. 



3. CALCULATING THE EIGENPROJECTION OF A BY 
ANY ANNIHILATING POLYNOMIAL FOR yl" 

The main result of this section is formulated as Theorem 1 stating that the eigenprojection of a 
matrix A can be found by direct calculation using any annihilating polynomial for any matrix A^, 
where u> v = indA. One can always take v as, u (see (1)) and the characteristic polynomial^ as 
the annihilating polynomial, but to reduce computations, it would be advantageous to make use of 
the information available about A for getting a more accurate upper bound for v and, perhaps, for 
obtaining an annihilating polynomial of a power lower than that of the characteristic polynomial. 
In this respect, the minimum polynomial for A'^ is the best choice, but its determination can be 
problematic itself. 

The proof of Theorem 1 will rely on the above Condition (f) which is necessary and sufficient for 
a matrix Z to be the eigenprojection of A. We note, however, that direct reduction to the definition 
of the eigenprojection or to other necessary and sufficient conditions differs insignificantly in labor 
consumption. 

Let, as before, u > v (for example, u = v oy u = n) and let ip{X) be an arbitrary nonzero 
polynomial which is annihilating for A^: ip{A^) = 0. Let 

<^(A) = A*(A«+PiA«-i + ...+;>,), (10) 

where t, g G {0, 1, . . .}, pi, . . . ,Pq G C, and pq ^ 0. 

We set 

h{X)=p-\X1+piX1-' + ...+Pq) (11) 

and 

Z = /i(A"). (12) 

^ The coefRcients of the characteristic polynomial (6) and the matrix coefEcients Ai of the adjoint matrix 
(7) can be calculated concurrently by Faddeev's algorithm [4,21]. 

It can be calculated using the Leverrier-Faddeev algorithm, the Krylov method, or other computational 
algorithms (see, for example, [4,21]). 
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Proposition 1. For any matrix Z determined by (12), the following statements hold: 

(1) ZC = CZ for any matrix C commuting with A; 

(2) Z is idempotent; 

(3) AZ is nilpotent; 

(4) if A is singular, then Z 0. 

Proof of Proposition 1. Item (1) follows from the fact that Z is & polynomial of A. 

(2) Let us consider the polynomial 

^o(A) = Xh{X) = p-U(A^ +piA«-i + . . . +pg). 

Let tp{\) be the minimum polynomial for A'^. Then, according to Theorem 3.3.1 of [25], ^^(A) is 
divisible by and, by virtue of Theorem 3.3.6 of [25] and Eq. (3), the lower power A in V'(A) is 
at most 1. Therefore, V'(A) divides v^o(A) and, consequently, V'o(A) is an annihilating polynomial 
for A^, that is, 

= ipo{A'^) = = A"Z = Z^". (13) 

Let us prove that 

Z^ = Z. 

Indeed, by using (12) and (13) we obtain 

Z^ = ZhiA"") = Z{A''g{A'') + I) = Z, 

where h{X) is represented as A g{\) + 1 and g{X) is a polynomial of A. 

(3) If n = 0, then Z = by virtue of (13); consequently, AZ is nilpotent with the nilpotency 
index 1. If n > 0, then by virtue of Items (1) and (2) and also Eq. (13) we obtain 

(AZ)'^ = A^Z"" = A^Z = 0, 

and AZ is nilpotent with the index of nilpotency at most u. 

(4) If A is singular, then according to Theorem 4.6 (c) of [11] the smallest power of A in the 
minimum polynomial V'(A) is equal to the index of A", that is, is strictly positive. Therefore, h{X) is 
not divisible by '0(A) and, consequently, h{X) is not an annihilating polynomial for A". Therefore, 

^ hiA") = Z. 



The following theorem results from Proposition 1 and Condition (f). 

Theorem 1. Let v = indA and u > v. Then, for any nonzero polynomial ^{X) which is 
annihilating for A^ , the matrix Z defined by (12) is the eigenprojection of A. 

The eigenprojection of A is related with the Drazin inverse A^ (see (g) in Sec. 2). If the 
eigenprojection is known, then the matrix A^ can be determined [20] by the formula^ 

A^ = {A + aZ)-\l-Z), a^O. (14) 
^ For a = 1, formula (14) was obtained in [15]. 
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Additionally, if ind^ = 1, then the group inverse A* = A^ is representable as 



A* = (A + Z)-^ -Z 

(see, for example, [6, 23]). Together with the formula A^ = A'^^^{A'^)* that can be easily checked 
and the fact that ind(^'') = 1, the above representation allows one to obtain readily another 
expression for A^: 

A^ = A''-\iA'' + Z)-^ - Z). 

Wc also note that the Moore-Penrose generalized inverses can also be expressed in terms of eigen- 
projections (see, for example, pp. 686-687 in [24]). 

4. DETERMINING THE COMPONENTS AND THE MINIMUM POLYNOMIAL OF A 

MATRIX 

Let \i, . . . ,Xs G C be all distinct eigenvalues of A. We denote by Vk the index of the eigenvalue 
Ajk , k = 1, . . . ,s, that is, the order of the greatest Jordan block corresponding to Afc, which is 
known to be equal to the power of the multiplier (A — A^) in the minimum polynomial of A. By 
the eigenprojection of A corresponding to A^ is meant the eigenprojection of the matrix A — X^I- 

Having determined Ai,...,As, one can use Theorem 1 to determine the corresponding eigen- 
projections. Now, we turn to the theory of matrix components [4, 5] according to which, for any 
function / : C — > C having finite derivatives f^^\\k) at the points Ai, . . . , Ag of the corresponding 
orders j = 0, . . . , f^, the function f{A) is naturally defined as follows: 

f{A) :=^^/(^)(A,)Z,„ 

k=l j=0 

where by the derivative of order is meant the value of the function. 

The matrices Zj^j are called the components of A. Here, the component Z^q is the eigenprojection 
of A corresponding to Afc, k = 1, . . . ,s, the rest of the components being obtained by successive 
multiplication of Z^q hy A — Xf^I along with some numerical coefficients: 

Zkj = {j\r\A-Xkiy Zko (15) 

(see Theorem 5.5.2 in [5]). 

All components are linearly independent (see, for example § 5.3 in [4]), thus, in particular, there 
are no zero components among them. At the same time, {A — Xk lY'' Z^o = for each k = 1, . . . , s, 
which follows from condition (b) of Sec. 2 and the fact that Z^q is the eigenprojection of ^ — A^ /, 
Vk being its index. Therefore, upon determining the component Zj^o by means of Theorem 1, the 
remaining components can be determined by multiplying Zi-q by ^ — Afc 7 and the corresponding 
numerical coefficients (see (15)). As soon as the next multiplication by A — X^I provides 0, all 
components corresponding to Xk are already determined and their mimbcr is Vk- 

In turn, knowledge of A^ and i>k allows one to construct the minimum polynomial of A as follows: 

m = f[{x-Xkr. 

k=l 

Thus, starting with the eigenvalues of A and some annihilating polynomials for 

(^-Afc)", u>mdiA-Xk), k = l,...,s, 
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one can further determine the eigenprojections corresponding to all eigenvalues of A, the compo- 
nents of A, the indices of the eigenvalues, and the minimum polynomial for A. 

The following problem arises when this method is realized numerically. The algorithm involves 
verification of whether the current multiplication of (A—X^ ly Z^q hy A—\j^I results in 0. Approx- 
imate calculations need an additional criterion for recognizing a matrix with "small" entries as the 
zero matrix. However, these problems arise with all algorithms for computing the rank of a matrix 
or the index of its eigenvalues because these values are unstable to arbitrarily small perturbations 
of the matrix entries. 

5. EIGENPROJECTION AND COMPONENTS OF A MATRIX 
WITH KNOWN EIGENVALUES 

One can readily check that if J {a) is a Jordan block with eigenvalue a 7^ 0, then ind {J{a) — al) = 
ind iJ^{a) — aPI), p = 2,3, . . . (see, for example, Problem 16 of Sec. 3.2 in [25]). This gives rise to 
the following lemma. 

Lemma 2. Let Aj ^ 6e the eigenvalue of A. Then, the index of Xi is equal to the index of Xf 
as the eigenvalue of A^, p = 2,3, 

The following expression for the eigenprojection and components of a matrix is useful if one 
knows the eigenvalues of A, rather than an annihilating polynomial of A. 

Proposition 2. Let A G C"^"; let v and Z be the index and eigenprojection of A, respectively; 
let Ai, . . . , As be its different eigenvalues, ui,. . . ,Us be their indices, and the integers u,ui, . . . ,Us 
be such that u> v and Ui > Vi, i = 1, . . . , s. Then, 

n iA--xpr^ 

Z = ^ yr- (16) 

n -Ar 

and 

n HA-XkiY^ -{\-XkTHf' 

Zkj = — J yr, ijr\A - Xkiy, (17) 

n -(A^-Afc)«o 

Wk J 

where Z^j is the order j component of A corresponding to X^, k = 1, . . . , s, j = 0, . . . — 1. 

For j = 0, Eq. (17) defines the eigenprojections of A corresponding to its eigenvalues. 
Proof of Proposition 2. We note that 

AC(A) = A n i^-Kr 

is an annihilating polynomial for A^. Indeed, this polynomial is divisible by the minimum polyno- 
mial for A^ because A" are the eigenvalues of A^, their indices do not exceed the numbers Ui by 
Lemma 2, and ind (yl") < 1 according to (3). To prove (16), it suffices now to take h{X) as the 
polynomial ^(A) divided by its absolute term 

p = ( n (-^D 
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and apply Theorem 1. 

By applying now (16) to the matrix A — X^I whose index does not exceed by definition and 
the distinct eigenvalues are Aj — A^, i = 1, . . . , s, we obtain (17) with j = 0. Now, the expressions 
(17) for the components of A of higher orders follow from (15). Proposition 2 generalizes formula 
(5.4.3) in [5] which refers to the case of ui = . . . = i/g = 1. 

6. CONCLUSIONS 

A number of necessary and sufficient conditions defining the eigenprojection (the principal idem- 
potent) of a square matrix were reviewed. A theorem which enables one to determine the eigen- 
projection of a matrix A and thereby the Drazin inverse A^ by means of any nonzero annihilating 
polynomial for any matrix A^ with u > ind^ was proved and used to determine the components 
of A and its minimum polynomial. 
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